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This article proposes a characterisation of the double-reed in quasi-static regimes. The non- 
linear relation between the pressure drop Ap in the double-reed and the volume flow crossing it q is 
measured for slow variations of these variables. The volume flow is determined from the pressure drop 
in a diaphragm replacing the instrument's bore. Measurements are compared to other experimental 
results on reed instrument exciters and to physical models, revealing that clarinet, oboe and bassoon 
quasi-static behavior relies on similar working principles. Differences in the experimental results are 
interpreted in terms of pressure recovery due to the conical diffuser role of the downstream part of 
double-reed mouthpieces (the staple). 

PACS numbers: 43.75.Pq 



I. INTRODUCTION 



Context 



The usual method for studying and simulating the be- 
havior of self-sustained instruments is to separate them 
in two functional parts that interact through a set of 
linked variables: the resonator, typically described by 
linear acoustics, and the exciter, a non-linear element in- 
dispensable to create and maintain the auto-oscillations 
of the instrument |l6j . 

Although this separation may be artificial because 
there may not be a well-defined boundary between the 
two systems, it is usually a simplified view that allows 
to describe the basic functioning principles of the instru- 
ment. In reed instruments, for instance, the resonator 
is assimilated to an air column inside the bore, and the 
exciter to the reed, which acts as a valve. 

In the resonator of wind instruments, the relation be- 
tween the acoustic variables, pressure (p) and volume 
flow (q) can be described by a linear approximation to 
the acoustic propagation which has no perceptive conse- 
quences in sound simulations ^f|. The oscillation arises 
from the coupling between the reed and the air column, 
which is mathematically established through the vari- 
ables p and q. 

The exciter is necessarily a nonlinear component, so 
that the continuous source of energy supplied by the pres- 
sure inside the musician's mouth can be transformed into 



an oscillating one The characterisation of the 

exciter thus requires the knowledge of the relations be- 
tween variables p and q at the reed output (the coupling 
region) . 

In principle this relation is non-instantaneous, because 
of inertial effects in the reed oscillation and the fluid dy- 
namics. A first approach to the characterisation of the 
exciter is thus to restrict the study to cases where these 
delayed dependencies (or, equivalently, time derivatives 
in the mathematical description of the exciter) can be 
neglected. 

The aim of this paper is to measure the relation be- 
tween the pressure and volume flow at the double reed 
output in a quasi-static case, that is, when the time varia- 
tions of p and q can be neglected, and to propose a model 
to explain the measured relation. 



B. Elementary reed model 

In quasi-static conditions, a simple model can be used 
to describe the reed behavior |2^| . The reed opening area 
(S) is controlled by the difference of pressure on both 
sides of the reed: the pressure inside the reed (p r ) and the 
pressure inside the mouth (p m )- In the simplest model, 
the relation between pressure and reed opening area is 
considered to be linear and related through a stiffness 
constant (k s ): 



(Ap) r = p m - p r = k s (S Q - S) 



(1) 
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In this formula, So is the reed opening area at rest, when 
the pressure is the same on both sides of the reed. In 
most instruments (such as clarinets, oboes or bassoons) 
the reed is said to be blown-closed (or inward- striking) 
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|l6j . because when the mouth pressure (p m ) is increased, 
the reed opening area decreases. 

The role of the reed is to control and modulate the 
volume flow (q) entering the instrument. The Bernoulli 
theorem applied between the mouth and the reed duct 
determines the velocity of the flow inside the reed (u r ) 
independently of the reed opening area: 



1 2 1 2 

Pm + ^pu m = Pr + -PU T 



(2) 



In this equation, p is the air density. Usually, the flow 
velocity u m is neglected inside the mouth, because of 
volume flow conservation: inside the mouth the flow is 
distributed along a much wider cross-section than inside 
the reed duct. 

The volume flow (q) is the integrated flow velocity (u r ) 
over a cross-section of the reed duct. For the sake of 
simplicity, the flow velocity is considered to be constant 
over the whole opening area, so that q — Su r . Using 
equation J5J), the flow is given by 



q = S\ 



'2(p m -p r ) 



P 



(3) 



Combining eq. and eq. it is possible to find the 
relation between the variables that establish the coupling 
with the resonator (p r and q): 



q 



PM ~ (Ap) r 2(Ap) r 

k s V P 



(4) 



The relation defined by equation Q is plotted in figure 
^ constituting what we will call the elementary model for 
the reed. 
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Equation Q can then be re-written in terms of these 
dimcnsionlcss quantities: 



'1 



(1 -P)P 



(7) 
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This formula shows that the shape of the non-linear 
characteristic curve of the elementary model is indepen- 
dent of the reed and blowing parameters, although the 
curve is scaled along the pressure p and volume flow q 
axis both by the stiffness k s and the beating pressure 
Pm = k s S . 



C. Generalisation to double-reeds 

For reed instruments, the quasi-static non-linear char- 
acteristic curve has been measured in a clarinet mouth- 
piece @, ^3] j and the elementary mathematical model 
described above can explain remarkably well the obtained 
curve almost until the reed beating pressure (pm). 

For double-reed instruments it was not verified that 
the same model can be applied. In fact, there are some 
geometrical differences in the flow path that can con- 
siderably change the theoretical relation of equation J7J). 
Local minima of the reed duct cross-section may cause 
the separation of the flow from the walls and an addi- 
tional loss of head of the flow j2^| , and in that case the 
characteristics curve would change from single-valued to 
multi- valued in a limited pressure range. This kind of 
change could have significant consequences on the reed 
oscillations. 

However, the non-linear characteristic relation was 
never measured before for double-reeds, justifying the 
work that is presented below. 



II. PRINCIPLES OF MEASUREMENT AND 
PRACTICAL ISSUES 

The characteristic curve requires the synchronised 
measurement of two quantities: the pressure drop across 
the reed (Ap) r and the induced volume flow q. 



FIG. 1: A typical non-linear characteristic curve for a reed of 
dimensions similar to an oboe reed, given by equation Q 



Volume flow measurements 



The static reed beating pressure pu = k s S [24[ 
is an alternative parameter to So, and can be used as a 
magnitude for proposing a dimcnsionlcss pressure: 



P = {Ap)r/PM 



(5) 



Similarly a magnitude can be found for q, leading to the 
definition of the dimcnsionlcss volume flow 



(G) 




One of the main difficulties in the measurement of the 
reed characteristics lies in the measurement of the vol- 
ume flow. There arc instruments which can accurately 
measure the flow velocity in an isolated point (LDA, hot- 
wire probes) or in a region of a plane (PIV), but it can 
be difficult to calculate the corresponding flow by inte- 
grating the velocity field. In fact it is difficult to do a 
sampling of a complete cross-section of the reed because 
a large number of points would have to be registered. 
Supposing that the flow is axisymmetric at the reed out- 
put (which is confirmed by experimental results Q), the 
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measurement along a diameter of the reed would be suf- 
ficient, but regions close to the wall are inaccessible. 

On the other hand, commercial flow-meters usually 
have the disadvantage of requiring a direct reading, which 
would have been unpractical for a complete characteristic 
measurement (large number of readings in a short time 
interval). 

An indirect way of measuring the flow was then pre- 
ferred to the above mentioned methods. It consists in 
introducing a flow resistance in series with the reed, for 
which the pressure can be accurately related to the flow 
running through it (see fig- - 

The diaphragm method, used successfully by S. 01- 
livier ^| to measure the non-linear characteristic of sin- 
gle reeds, is based on this principle. The resistance is 
simply a perforated metal disk which covers the reed out- 
put. 



(A,,). 




FIG. 2: Use of a diaphragm to measure flow and pressure 
difference in the reed. Numbered rectangles correspond to 
the pressure probes used in the measurement. 

For such a resistance, and assuming laminar, viscous 
less flow, the pressure drop (Ap) d = p r —p a tm across the 
diap hragm can be approximated by the Bernoulli law 

m 



(Ap) 



Pr Patm 



1 5 
2'U 



(8) 



where q is the flow crossing the diaphragm, Sd the cross 
section of the hole, and p the density of air. In our experi- 
ment, pressure p a tm is the pressure downstream of the di- 
aphragm (usually the atmospheric pressure, because the 
flow opens directly into free air). The volume flow q is 
then determined using a single pressure measurement p r . 



B. Practical issues and solutions 

1. Issues 

The realisation of the characteristic measurement ex- 
periments encountered two main problems: 

a. Diaphragm reduces the range of (Ap) r for which 
the measurement is possible The addition of a resistance 
to the air flow circuit of the reed changes the overall non- 
linear characteristic of the reed plus diaphragm system 
(corresponding to (Ap) s in fig. |3 and to the dashed line 



in fig. |2J). The solid line plots the flow against (Ap) , 
the pressure drop needed to plot the non-linear charac- 
teristics. When the resistance is increased, the maximum 
value of the system's characteristic is displaced towards 
higher pressures , whereas the static beating pressure 
(pm) value does not change (because when the reed closes 
there is no flow and the pressure drop in the diaphragm 
(( A P)d) is zero)- 



q m 




FIG. 3: Comparison of the theoretical reed characteristics 
(solid line) with the a model of the overall characteristics of 
the reed associated to a diaphragm (dashed) — mathematical 
models, based on the Bernoulli theorem 
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Therefore, if the diaphragm is too small (i. e., the 
resistance is too high), part of the decreasing region 
(B'C) of the system's characteristics becomes vertical, 
or even multivalued, so that there is a quick transition 
between two distant flow values, preventing the measure- 
ment of this part of the characteristic curve [13 (see fig- 
ure |3j). A critical diaphragm size {Sd, C rit = 0.585o) can 
be found below which the characteristic curve becomes 
multi- valued (see appendix |A"|) . 

b. Reed auto- oscillations Auto-oscillations have to 
be prevented here to keep consistent with the quasi-static 
measurement (slow variations of pressure and flow) . This 
proved to be difficult to achieve in practice. In fact, auto- 
oscillations become possible when the reed ceases to act 
as a passive resistance (a positive J^, which absorbs en- 
ergy from the standing wave inside the reed channel) to 
become an active supply of energy (^ < 0). All real 
acoustic resonators are slightly resistive (the input ad- 
mittance Yi n has a positive real part). This can com- 
pensate in part the negative resistance of the reed in its 
active region, but only below a threshold pressure, where 
the slope of the characteristic curve is smaller than the 
real part of Yi„ for the resonator ^3] ■ 

One way to avoid auto-oscillations is thus to increase 
the real part of Yi n , that is the acoustic resistance of the 
resonator. It is known that a an orifice in an a acoustical 
duct with a steady flow works as an anechoic termina- 
tion works as an acoustic resistance 01 > so that if the 
diaphragm used to measure the flow (see sect. Ill All is 
correctly dimensioned, the acoustic admittance seen by 
the reed Y in can become sufficiently resistive to avoid 
oscillations. 
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2. Solutions proposed to adress these issues 

a. Size of the diaphragm The volume flow is de- 
termined from the pressure drop across the diaphragm 
placed downstream of the reed. In practice there's a 
tradeoff that determines the ideal size of the diaphragm. 
If it is too wide, the pressure drop is too small to be mea- 
sured accurately, and reed oscillations likely to occur. If 
the diaphragm is too small, the system-wide character- 
istic can become too steep, making part of the ((Ap) r ) 
range inaccessible. 

The ideal diaphragm cross-section is then found em- 
pirically, by trying out several resistance values until one 
complete measurement can be done without oscillations 
or sudden closings of the reed. The optimal diaphragm 
diameter is sought using a medical flow regulator with 
continuously adjustable cross-section as a replacement for 
the diaphragm. 

b. Finer control of the mouth pressure p rn During 
the attempts to find an optimal diaphragm, it was found 
that sudden closures were correlated to sudden increases 
in the mouth pressure. A part of the problem is that 
the mouth pressure depends both on the reducer setting 
[or configuration] and on the downstream resistance. By 
introducing a leak upstream of the experimental appa- 
ratus (thus not altering the experiment) , it is possible 
to improve measurements in the decreasing region of the 
characteristic (BC), at least when the system-wide char- 
acteristic is not multi- valued (see section Hi B 1 a|) . 

c. Increase the reed mass One other way to reduce 
the oscillations is thus to prevent the appearance of in- 
stabilities, or to reduce their effects. An increase in the 
reed damping would certainly be a good method to avoid 
oscillations, because it cancels out the active role of the 
reed (which can be seen as a negative damping) 

It is difficult to increase the damping of the reed with- 
out altering its opening or stiffness properties. The sim- 
plest way found to prevent reed oscillations was thus an 
increase in the reed mass. 

This mass increase was implemented by attaching 
small masses of Blu- Tack pel] to one or both blades of the 
reed (fig-EJ. During measurements on previously soaked 
reeds it was difficult to keep the masses attached to the 
reed, so that an additional portion of Blu-Tack is used 
to connect the to masses together, wrapping around the 
reed. A comparison of the results using different masses 
showed that their effect on the quasi-static characteristics 
can be neglected. 

C. Experimental set-up and calibrations 

The experimental device is shown in figure An ar- 
tificial mouth was used as a blowing mechanism and 
support for the reed. The window in front of the reed al- 
lows the capture of frontal pictures of the reed opening. 

Artificial lips, allowing to adjust the initial opening 
area of the reed were not used here, fearing that they 




FIG. 4: Front view of the reed with attached masses, at left 
in dry conditions, at right in soaked conditions (to prevent 
the masses from slipping). 



would modify some of the clastic properties of the reed, 
yet differently from what happens with real lips. 




FIG. 5: Device used for characteristics measurements. 

As stated before, the plot of the characteristic curve re- 
quires two coordinated measurements: the pressure dif- 
ference (Ap) r across the reed and the induced volume 
flow q, determined from the pressure drop (Ap)d across 
a calibrated diaphragm (sect. Ill AT) . 

In practice thus, the experiment requires two pressure 
measurements p m and p r , as shown in figure [SJ 

1. Pressure measurements 

The pressure is measured in the mouth and in the reed 
using Honeywell SCX series, silicon-membrane differen- 
tial pressure sensors whose range is from —50 to 50 kPa. 

These sensors are not mounted directly on the mea- 
surement points, but one of the terminals in each sen- 
sor is connected to the measurement point using a short 
flexible tube (about 20 cm in length). Therefore, one 
tube opens in the inside wall of the artificial mouth, 4 
cm upstream from the reed, and the other tube crosses 
the rubber socket attaching the diaphragm to the reed 
output. The use of these tubes does not influence the 
measured pressures as long as their variations are slow. 

The signal from these sensors is amplified before en- 
tering the digital acquisition card. The gain is adjusted 
for each type of reed. The system consisting of the sen- 
sor connected to the amplifier is calibrated as a whole 
in order to find the voltage at the amplifier output cor- 
responding to each pressure difference in the probe ter- 
minals: the stable pressure drop applied to the probe is 
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also measured using a digital manometer connected to 
the same volumes, and compared to the probe tension 
read using a digital voltmeter. Voltage is found to vary 
linearly with the applied pressure within the measuring 
range of the sensor. 



2. Diaphragm calibration 

The curve relating volume flow q to the pressure dif- 
ference through diaphragms (A p )d can be approximated 
by the Bernoulli theorem. In fact, diaphragms are con- 
structed so as to minimize friction effects (by reducing 
the length of the diaphragm channel) and jet contraction 
- the upstream edges are smoothed by chamfering at 
45° (fig. EJ). The chamfer height (c) is approximately 0.5 
mm. The diaphragm channel is 3 mm long (L). 




to turbulence generated for high Reynolds Numbers. In 
figureElthe dashed line corresponding to the critical value 
of the Reynolds number (Re c = ^7 = 2000) is shown. It 
is calculated using the following formulas for u (the aver- 
age flow velocity in the diaphragm) and d (the diaphragm 
diameter) : 



7T U 



1/2 



2(Ap) d 



1/2 



so that the constant Reynolds relation is given by: 



Q 1/2 Ap 1/4 



Re c v ( ~ 
V2 



1/2 



1/4 



(9) 
(10) 

(11) 



where the right-hand side should be a constant based on 
the diaphragm geometry. 

Since a suitable model was not found for the data dis- 
played in figure[7| we chose to interpolate the experimen- 
tal calibrations in order to find the flow corresponding to 
each pressure drop in the diaphragm. Linear interpola- 
tion was used in the (p, q 2 ) space. 



FIG. 6: Detail of the diaphragm dimensions. 



3. Typical run 
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FIG. 7: Calibration of diaphragms used in characteristic mea- 
surements (dots are experimental data and lines are Bernoulli 
predictions using the measured diaphragm diameters). The 
dashed black line represents the pressure/flow relation corre- 
sponding to the expected transition between laminar and tur- 
bulent flows (Re c = 2000), parameterised by the diaphragm 
diameter d. 

Nevertheless, this ideal characteristic was checked for 
each diaphragm (see fig [7J . It was found that the ef- 
fective cross-section is slightly smaller than the actual 
cross-section (about 10%), which is probably due to some 
Vena Contracta effect in the entrance of the diaphragm. 
Moreover, above a given pressure drop the volume flow 
increase is lower than what is predicted by Bernoulli's 
theorem (corresponding to a lower exponent than 1/2 
predicted by Bernoulli). This difference is probably due 



In a typical run, the mouth pressure p m is equilibrated 
with the atmospheric pressure in the room in the begin- 
ning of the experiment. Both p m and p r are recorded 
in the computer through a digital acquisition device at a 
sampling rate of 4000 Hz. p m is increased until slightly 
above the pressure at which the reed closes, left for some 
seconds above this value and then decreased back to 
the atmospheric pressure. The whole procedure lasts for 
about 3 minutes, and is depicted in figure |SJ 
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FIG. 8: Time variation of the mouth pressure (p m )and the 
pressure inside the reed (p r ) during a successful characteristic 
measurement. 
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D. Double- reeds used in this study and operating 
conditions 



III. RESULTS AND DISCUSSION 



A. Typical pressure vs flow characteristics 



Among the great variety of double-reeds that are used 
in music, we chose as a first target for these measure- 
ments a natural cane oboe-reed fabricated using standard 
procedures (by Glotin), and sold to the oboist (usually 
a beginner oboist) as a final product (i.e., ready to be 
played). 

The choice of a ready-to-use cane reed was mainly re- 
tained because it can be considered as an average reed. 
This avoids considering a particular scraping technique 
among many used by musicians and reed-makers. Of 
course, this does not greatly facilitate the task of the reed 
measurement, because natural reeds are very sensitive to 
environment conditions, age or time of usage. 

Other reeds were also tested, as a term of comparison 
with the natural reeds used in most of the experiments. 
However, none of these reeds was produced by a pro- 
fessional oboist or reed maker, although it would be an 
interesting project to investigate the variations in reeds 
produced by different professionals. 

To conclude, the results presented in the next section 
may depend to a certain extent on the reed chosen for 
the experiments, and a larger sample of reeds embracing 
the big diversity of scraping techniques needs to be tested 
before claiming for the generality of the results that will 
be presented. 

Another remark has to be made on the conditions dur- 
ing the experiments. The kind of reeds used in most 
experiments are always blown with highly moisturised 
air. In fact, in real life, reeds are often soaked before 
they are used, and constantly maintained wet by saliva 
and water vapour condensation. These conditions were 
sought throughout most of the experiments, although 
their sensitivity to environmental conditions was also in- 
vestigated. For instance, the added masses were found to 
have no practical influence on the nonlinear characteris- 
tics, whereas the humidity increases the hysteresis in the 
complete measurement cycle (increasing followed by de- 
creasing pressures), while reducing the reed opening at 
rest 0. 

In our measurements, humidification is achieved by let- 
ting the air flow through a plastic bottle half-filled with 
hot water at 40° (see fig. El, recovering it from the top. 
Air arriving in the artificial mouth has a lower tempera- 
ture, because its temperature is approximately 10° when 
entering the bottle. This causes the temperature and hu- 
midity to decrease gradually along the experiments. Fu- 
ture measurements should include a thermostat for the 
water temperature in order to ensure stable humidifica- 
tion 



Using the formula of eq. (JSJ, and the calibrations car- 
ried out for the diaphragm used in the measurement, 
the flow is determined from the pressure inside the reed 
(p r ). The pressure drop in the reed corresponds to 
the difference between the mouth and reed pressures 
((Ap) r — p m — p r ). Volume flow (q) is then plotted 
against the pressure difference ((Ap) r ), yielding a curve 
shown in figure 
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FIG. 9: A typical result for the measurement of the volume 
flow vs pressure characteristic of a natural cane oboe reed. 

In this figure, the flow is seen to increase until a certain 
maximum value (at about 6 kPa). When the pressure is 
increased further, flow decreases due to the closing of the 
reed. Instead of completely vanishing for (Ap) r = pm, as 
predicted by the elementary model shown in section IT*Bl 
the flow first stabilises at a certain minimum value and 
then increases when the pressure is increased further., 
indicating that it is very hard to completely close the 
reed. 

The flow remaining after the two blades are in contact 
suggests that despite the closed apparency of the double 
reed, some narrow channels remaining between the two 
blades are impossible to close, behaving like rigid capil- 
lary ducts, which is corroborated by the slight increase 
in the residual volume flow for high pressures. 

If the residual reed opening is distributed over several 
channels, the flow in each of these channels would be con- 
trolled by viscosity rather then inertia, so that the volume 
flow q should be proportional to (Ap) r rather than to the 
1/2 power of it as in the case of a Bernoulli flow. Nev- 
ertheless, observing a logarithmic plot of the non-linear 
characteristics ffig. HOf) shows a 1/2 power dependence of 
the residual flow. This may mean that the reed material 
is not yet stable at this point. 

When reducing the pressure back to zero, the reed fol- 
lows a different path in the pressure/flow space than the 
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Pressure flow characteristic 
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FIG. 10: Logarithmic plot of the characteristic curve of figure 
[§|to show the 1/2 power dependence when the reed is almost 
shut. 



path for increasing pressures. This hysteresis is due to 
memory effects of the reed material which have been in- 
vestigated experimentally for single-reed [Tol ] and double- 
reed instruments ||. 



B. Comparison with other instruments 

1. Bassoon 

Since oboes are not the only doublc-recd instruments, 
it is interesting to compare the non-linear characteristic 
curves from different instruments. The bassoon is also 
played using a double reed, but its dimensions are differ- 
ent: its opening area at rest is typically 7 mm 2 (against 
around 2 mm 2 for the oboe) and the cross section profile 
varies slightly from oboe reeds. 

Figure 1111 compares two characteristic curves for nat- 
ural cane oboe and bassoon reeds. Both were measured 
under similar experimental conditions, as far as possible. 
The reed was introduced dry in the artificial mouth, but 
the supplied air is moisturized at nearly 100% humidity, 
and masses were added to both reeds to prevent auto- 
oscillations. The diaphragms used in each measurement 
are different however, and this is because the opening 
area of the reed at rest is much larger in the case of the 
bassoon, so that a smaller resistance (a larger diaphragm) 
is needed to avoid that the reed closes suddenly in the de- 
creasing side of the characteristic curve (see section III CI 
and eq. IA10|) . This should not have any consequences in 
the measured characteristic curve. 

In the q axis, the bassoon reed reaches higher values, 
and this is probably a consequence of its larger opening 
area at rest, although the surface stiffness is likely to 
change as well from the oboe to the bassoon reed. In the 
p axis, the bassoon reed extends over a smaller range of 
pressures so that the reed beating pressure is about 17 
kPa in the case of the bassoon reed whereas it is near 33 



FIG. 11: Comparison of the characteristic curves of differ- 
ent reed exciters for different instruments. Clarinet data 
was obtained by S. Ollivier and J. -P. Dalmont [l(J for a 
PlastiCover® reed. Oboe and bassoon reeds used moisturised 
air. 



kPa for the oboe reed. 

Apart from these scaling considerations, the shape of 
the curves are similar and this can be better observed 
if flow and pressure are normalized using the maximum 
flow point of each curve (figure EJ- 
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FIG. 12: Normalized data from figure fTTI 



2. Clarinet 

The excitation mechanism of clarinets and saxophones 
share the same principle of functioning with double reeds. 
However, there are several geometric and mechanical dif- 
ferences between single reeds and double reeds. For 
instance, flow in a clarinet mouthpiece encounters an 
abrupt expansion after the first 2 or 3 millimeters of the 
channel between the reed and the rigid mouthpiece, and 
the single reed is subject to fewer mechanical constraints 
than any double reed. These differences suggest that 
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the characteristic curve of single-reed instruments might 
present some qualitative differences with respect to the 
double reed [2Cj . 

The non-linear characteristic curve of clarinet mouth- 
pieces displayed in figures ITT1 and IT2l was measured by S. 
Ollivier an J. -P. Dalmont [llj using similar methods as 
the ones we used for the double reed. A comparison be- 
tween the curves for both kinds of exciters (in figure ITT|) 
shows that the overall behavior of the excitation mech- 
anism is similar in both cases. Similarly to when com- 
paring oboe to bassoon reeds, the scalings of the charac- 
teristic curves of singlc-rceds arc different from those of 
oboe reeds, although closer to those of the bassoon. This 
is probably a question of the dimensions of the opening 
area. 

A different issue is the relation between reference pres- 
sure values in the curve (shown in the adimensionalizcd 
representation of figure ll2|) . As predicted by the elemen- 
tary model described in section lFBl in the single reed the 
pressure at maximum flow is about 1/3 of the beating 
pressure of the reed, whereas in double-reed measure- 
ments, the relation seems to be closer to 1/4. This devi- 
ation from the model is shown in section Hvl to be linked 
with the diffuser effect of the conical staple in double- 
reeds. 

Figure El also shows that in the clarinet mouthpiece 
used by S. Ollivier and J. P. Dalmont the hysteresis is rel- 
atively less important than in both kinds of double-reeds. 
In fact, whereas the measurements for double-reeds were 
performed in wet conditions, the PlastiCover® reed used 
for the clarinet was especially chosen because of its 
smaller sensitivity to environment conditions. 
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pressure difference (kPa) 



50 



FIG. 13: Comparison of the experimental non-linear char- 
acteristics curve with the elementary model shown in figure 
Fwo models are fitted, for increasing (j>m = 35 kPa, 
k s = 1.04 x 10 10 kg m~ 3 s -2 ) and decreasing (p M = 27 kPa, 



k s = 8.86 x 10 a kg m" 



mouth pressures. 



reaches its maximum value. In fact, it can be easily 
shown that for the elementary model this value is 1/3 
of the reed beating pressure pm, which is also verified in 
the clarinet (sect. IHIB 2|) . In the measured curves how- 
ever, this value is usually situated between 1/ApM and 

Nevertheless, the shapes of the curves are qualitatively 
similar to the theoretical ones. 



B. Conical diffuser 



IV. ANALYSIS 
A. Comparison with the elementary model 

The measured non-linear characteristic curve of figure 
El can be compared to the model described in section llBl 
In this model, two parameters (k s and So) control the 
scaling of the curve along the p and q axis. They are 
used to adjust two key points in the theoretical curve to 
the experimental one: the reed beating pressure pm and 
the maximum volume flow q max . 

Once g m ax is determined through a direct reading, the 
stiffness k s is calculated using the following relation: 

ks = q m l x 1 p- 1 ' 2 (12) 

This allows to adjust a theoretical characteristic curve 
(corresponding to the elementary model of equation Q) 
to each of the branches of the measured characteristic 
curve, for increasing and decreasing pressures ffig. I13|) . 

When compared to the elementary model of section 
II Bl the characteristic curve associated to double-reeds 
shows a deviation of the pressure at which the flow 



The former observations about the displacement of the 
maximum value can be analysed in terms of the pres- 
sure recoveries due to flow decelerations inside the reed 
duct. Variations in the flow velocity are induced by the 
increasing cross-section of the reed towards the reed out- 
put ffig. ll4|) . This can be understood simply by consider- 
ing energy and mass conservation between two different 
sections of the reed: 

Pin + \p {-§-] = Pout + \p ( 7T- ) (13) 

* \^in J £ \^out / 

where q is the total volume flow that can be calculated 
either at the input or the output of the conical diffuser 
by integrating the flow velocity over the cross-section S in 
or S ou t respectively. 

In practice however, energy is not expected to be com- 
pletely conserved along the flow because of its turbulent 
nature. In fact, for instance at the reed output (diameter 
d), the Reynolds number of the flow (Re = ^ = 
can be estimated using data from figurc|^|to reach a max- 
imum value of 5000. Given that this number is inversely 
proportional to the diameter of the duct d, the Reynolds 
increases upstream, inside the reed duct, so that the flow 
is expected to be turbulent also for lower volume flows. 
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m 




FIG. 14: Cross-section profiles (axis and area) of an oboe 
reed, measured on a mold of the reed channel, and indexes 
used in section llV Bl mouth, constriction and reed output. 



For turbulent flows, no theoretical model can be ap- 
plied to calculate the pressure recovery due to the taper- 
ing of the reed duct. However, phenomenological models 
are available in engineering literature, where similar duct 
geometries are known as "conical diffusers" . Unlike in 
clarinet mouthpieces, where the sudden expansion of the 
profile is likely to cause a turbulent mixing without pres- 
sure recovery [17^ , this effect must be considered in coni- 
cal diffusers. The pressure recovery is usually quantified 
in terms of a recovery coefficient Cp stating the relation 
between the pressure difference between both ends of the 
diffuser and the ideal pressure recovery which would be 
achieved if the flow was stopped without losses: 



Pout - Ph 



(14) 



Cp values range from (no recovery) to 1 (complete re- 
covery, never achieved in practice). Distributed losses 
due to laminar viscosity along the reed arc neglected. 

According to equation (|14f> . pressure recovery is pro- 
portional to the square of the flow velocity at the en- 
trance of the conical diffuser, and consequently to the 
squared volume flow inside the reed. This is coherent 
with the shifting of the double-reed non-linear character- 
istics curve for high volume flows, as observed in figure 

EH 



1. Reed model with pressure recovery 

In order to take into account the pressure recovery be- 
fore the reed output, the flow is divided into two sections, 
the upstream, until the constriction at 28 mm (index c 
in fig. I14|) and the conical diffuser part from the constric- 
tion until the reed output. In the upstream section, no 
pressure recovery is considered, so that the flow velocity 
can be calculated using the pressure difference between 
the mouth and this point using a Bernoulli model, as in 
equation (j2J), but replacing p r with p c : 



q = S\ 



<2(p m -p c ) 



P 



(15) 



Similarly, the reed opening is calculated using the same 
pressure difference: 



(Ap) c = Pm - Pc = ks(S - S) 



(16) 



The total pressure difference used to plot the charac- 
teristic curve, however, is different, because the recovered 
pressure has to be added to (Ap) c : 

Pm ~Pr = (Ap) c - Cp-p ( -J- ) (17) 



where S c is the reed duct cross-section at the diffuser 
input i. e., at the constriction, which is found from figure 
d = 4 x 10" 6 m 2 . 

Using these equations, the modified model can be fit- 
ted to the experimental data. In figure 1151 the same 
parameters k s and So were used as in figure ITS1 leaving 
only C p as a free parameter for the fitting. Figure [TBI was 
obtained for a value of C„ = 0.8. 



0.35 




10 20 30 

pressure difference (kPa) 



FIG. 15: Comparison of the experimental non-linear charac- 
teristics curve with a reed model with pressure recovery in 
the final part of the duct. Fitted k s and So are the same as 
in figure Upl and different for increasing and decreasing mouth 
pressures. In both cases the value C p = 0.8 was used. 

This value can be compared to typical values of pres- 
sure recovery coefficients found in industrial machines, 
for example in Q. 
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In engineering literature, Cp is found to depend mostly 
on the ratio between output and input cross-sections 
(AR = S out /Sin) and the diffuser length to initial diame- 
ter ratio (L/d[ n ) |2l|. The tapering angle 6 influences the 
growth of the boundary layers, so that above a critical 
angle (6 — 8°) the flow is known to detach from the dif- 
fuser walls considerably lowering the recovered pressure. 
An in-depth study of turbulent flow in conical diffuscrs 
can be found in the literature usually for diffusers 
with much larger dimensions than the ones found in the 
double reed. 

The geometry of the conical diffuser studied in can 
be compared to the one studied in our work: the cross- 
section is circular and the tapering angle 6 — 3.94° is not 
very far from the tapering angle of the reed staple 6 = 
5.2° (in particular, both are situated in regions of similar 
flow regimes |l8j as a function of the already mentioned 
AR and L/d- m )- Reynolds numbers of his flows (Re = 
6.9 x 10 4 ) are also close to the maximum ones found at 
the staple input (Re ~ 10 4 ). 

The length to input diameter ratio of the reed staple 
L/d = 20 is bigger than that found in 0, however the 
pressure recovery coefficient can be extrapolated from his 
data to find the value C P ~ 0.8 (fig. 2 in 0). 



are quasi-static). Some clues indicate that this could also 
be put into question: 

• First of all, experimental observations of the flow 
at the output of the staple (through hot-wire mea- 
surements 0) revealed significant differences in the 
flow patterns when considering static and auto- 
oscillating reeds. 

• Moreover, non-dimensional analysis reveal Strouhal 
numbers much larger than for simple reed instru- 
ments (|H, 0) 
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APPENDIX A: CALCULATION OF THE 
MINIMUM DIAPHRAGM CROSS-SECTION 



CONCLUSION 



The total pressure drop in the reed-diaphragm system 
(fig. 13 is: 



The quasi-static non-linear characteristics was mea- 
sured for double reeds using a similar device as for single- 
reed mouthpieces [HI- The obtained curves are close to 
the ones found for single-reeds, and in particular no evi- 
dence of multivalued flows for a same pressure was found, 
as was suggested by some theoretical considerations |22| , 

M- 

Double-reed characteristic curves present substantial 
quantitative differences for high volume flows when com- 
pared to elementary models for the reed. These differ- 
ences can be explained using a model of pressure recovery 
in the conical staple, proportional to the input flow ve- 
locity. 

It is worth noting that direct application of this mea- 
surement to the modelling of the complete oboe is not 
that obvious. Indeed, the assumption that the mouth- 
piece as a whole (reed plus staple) can be modelled as 
a non-linear element with characteristics given by the 
above experiment would be valid if the size of the mouth- 
piece was negligible with respect to a typical wavelength. 
Given the 7 cm of the mouthpiece and the 50 cm of a typ- 
ical wavelength, this is questionable. Should the staple 
be considered as part of the resonator? In that case, the 
separation between the exciter and the resonator would 
reveal a longuer resonator and an exciter without pres- 
sure recovery. This could be investigated through nu- 
merical simulations by introducing the pressure recovery 
coefficient (Cp) as a free parameter. 

Moreover, the underlying assuption of the models is 
that non stationnary effects are negligible (all flow models 



(A P ) S = (A P ) r + (Ap) d 



(Al) 



The system's characteristics becomes multivalued 
when there is at least one point on the curve where the 
slope is infinite: 

A (Ap)s = | (Ap)r + ^ (Ap)d = (A2) 

Because of simplicity, the derivatives in equation (|A1|) 
are replaced by their inverse: 



yielding 



S /2(Ap), 



p V p 



Solving for Sd- 



-i/2 



1 /2(Ap) r 



k s \ p 



1/2N 



P /2(Ap) d 



Sd \ p 



1/2 



(A4) 
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Simplifying 

From equations (J2J and ©, we can find 



(A6) 



It is clear that the right-hand side of this equation 
must be positive. Moreover, it is a parabolic function of 



(AS) 



(Ap) r 



with its concavity facing downwards. 



The maximum value of S^(S): 



(f ) 



and equation (|A6() can be written 
S d 



_ S 2 S , . , 

3d i^s >Jd 



Now we can replace S = Sq — u to find 



5S = (-5b + 3^)(5b-^) 



(A7) 



(A8) 



(A9) 



max 



[S 2 d (S)) 



<? 2 
a 



(A10) 



is thus the value for which there is only a single point 
where the characteristic curve has an infinite slope. 



We thus conclude that Sd 



V3 



0.585*0 is the min- 



imum value of the diaphragm cross-section that should 
be used for flow measurements. 
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